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ABSTRACT 

In this study we present three-dimensional radiative cooling hydrodynamical simu- 
lations of galactic winds generated particularly in M82-like starburst galaxies. We have 
considered intermittent winds induced by SNe explosions within super star clusters 
randomly distributed (in space and time) in the central region of the galaxy (within 
a radius R = 150 pc) and were able to reproduce the observed M82 wind conditions 
with its complex morphological outflow structure. We have found that the environ- 
mental conditions in the disk in nearly recent past are crucial to determine whether 
the wind will develop a large scale rich filamentary structure, as in M82 wind, or not. 
If a sufficiently large number of super stellar clusters is build up in a starburst mainly 
over a period of a few million years, then the simulations reproduce the multi-phase 
gas observed in M82-like winds, i.e., with filaments of sizes about 20 to 300 pc, veloc- 
ities of ~ 200 - 500 km/s, densities in the range of 10 _1 - 10 cm~ 3 , embedded in a 
hot low density gas with a density smaller than 10~ 2 cm -3 and a velocity of ~ 2000 
km s _1 . Otherwise, a " superbubble-like" wind develops, with very poor or no cold 
filamentary structures. Also, the numerical evolution of the SN ejecta have allowed us 
to obtain the abundance distribution over the first ~ 3 kpc extension of the wind and 
we have found that the SNe explosions change significantly the metallicity only of the 
hot, low-density wind component for which we obtained abundances ~ 5 to 10 Z@ in 
fair consistency with the observations. Moreover, we have found that the SN-driven 
wind transports to outside the disk large amounts of energy, momentum and gas, but 
the more massive high-density component reaches only intermediate altitudes smaller 
than 1.5 kpc. Therefore, no significant amounts of gas mass are lost to the IGM and 
the mass evolution of the galaxy is not much affected by the starburst events occurring 
in the nuclear region. 
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1 INTRODUCTION 

The evolution of the interstellar medium (ISM) at galactic 
scales depends on the stellar feedback at the larger scales, 
and vice-versa. In fact, massive stars exert a profound influ- 
ence over the baryonic component of the galaxy. At the same 
time, the collective effect of supernova (SN) explosions, stel- 
lar winds and ra diation pressure may drive galactic winds 



lation pressure may driv e galactic winds 
Mathews fc Baked Tl97ll . IHeckman et al.1 



(GW s; 

Il99d . IVeilleux et alj 120051). GWs i nfluence the multi-phase 
conditions of the ISM l|McKeei 19951 ) as well as the dynamica 
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and chemical evolution of the galaxy and the thermodynam- 
ics and enrichment of the intergalactic medium (IGM; see, 
eg. 



Heckman et~alT l990. Shapiro et alJll994lAguirre et all 



12001 ). They are particularly important for starburst (SB) 
galaxies whose intrinsically large star formation rate (SFR) 
provides the perfect environment for gas outflow to develop. 

There has been extensive investigation in the literature 
on starburst GWs (see below). Our aim here is to extend 
these studies exploring the evolution of these winds in the 
surrounds of the galaxy. We will address questions which are 
not yet clearly understood such as: (i) how does the feedback 
of the stars of the nuclear region shape the wind? (ii) how 
much mass is thermally driven? (iii) how are the different 
phases of the gas connected? (iv) what is the fraction of cold 
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structures built up in the winds? (v) and how does the GW 
impact the chemical evolution of the galaxy? In order to try 
to answer to these questions, we have adopted the galaxy 
M82 as a reference to obtain appropriate initial conditions, 
although the results here presented may be applied to others 
GWs and outflows in general. 

M82 is a starburst galaxy, located at about 3.9 Mpc 
away, with a disk of ~ 5 kpc radius and mass ~ 10 10 
Mq. Although this galaxy has been classified as an Ir- 
regular type II, s e veral observations jMayva et alJ 120051 : 
iBarker et all l20Qg| ; IStrickland fc Heckmanl 120091 ') revealed 
the presence of a bar (~ 1 kpc) and spiral arms (observed 
in near infrared). The main feature of this galaxy is a 
GW with filamentary structures extending up to ~ 3 kpc 
above the disk and an Ha cap, i.e. a large nebular sur- 
rounding feature prob ably associated with the wind, at a 
height of about 10 kpc (fDevine fe B ally 199 91; iLehnert et al.l 
ll999l : lEngelbracht et alhood : IStrickland fc Heckmarj2009h . 
Many studies indicate that this galaxy has experienced a 
tidal interaction with the spiral galaxy M81, resulting a large 
amount of gas that is being piped into M82's core over the 
last ~ 200 Myr. The most recent and important interaction 
has happened ~ 2 - 5 x 10 8 yr ago, causing an intense star- 
burst which lasted about 50 Myr with a SFR ~ 10 M Q /yr. 
Then two more subsequent starbursts occurred, and the last 
one (about 4-6 Myr ago) may have for med at least one of the 
observed largest central star clusters ((Barker et al.ll2008h . 

The SB region in the core of M82 has a diameter of ~ 1 
kpc and is optically defined by por tions with high surface 
bright ness labelled A, C, D and E bv lO'Connell fc Mangand 
l|l978h . These regions correspond to sources of X-rays, in- 
frared and radio emission. The bipolar GW seems to be 
concentrated within A and C regions and is possibly being 
driven mostly by SNe energy injections. The SFR leads to a 
current supernovae rate ^ 0.1 year -1 . The central region has 
a mass of ~ 7 x 10 8 Mq, pressure ^ 10 -9 dyn cm -2 , tem- 
erature ~ 10 4 - 10 6 K and densities ~ 10 3 - 10 4 cm" 3 
1977; lO'Connell fc Mangand 1 19781 ; ISchaaf et al.1 



perature 
JCottrelll 



ll989l : IWestmoauette et alJl2007l .l 009al ) . The SNe inject mass 
and energy at rates of ~ 1 Mp yr -1 and ~ 2 x 10 42 erg/s, 
respectively (|Strickland fc Heckman|[2009l ) . Recent observa- 
tions revealed that within a radius of 500 pc there are about 
200 super stellar clusters (SSCs) with an average size of 
~ 5.7 pc and m ass (of stars) between 10 4 and 10 6 Mq 
ijMelo et al.ll2005l ). The projected separation between each 
of the young SSCs is smaller than 30 pc, with a minimum 
value of 5 pc and a mean separation of ~ 12 pc. The den- 
sity of young SSCs is very high in the M82 st a rburst galaxy, 
with a global value inferred by iMelo et al.1 l^OPd ) of 620 
kpc -2 , about 6 times larger than that observed in the nu- 
clear regions of ot hers SB galaxies (see, e.g., NGC 253, 
IWatson et aflll99d). The 20 brightest clusters have ages ~ 5 
to 10 x 10 6 yr (|McCradv et al.ll2003l ; IStrickland fc Heckmad 
l2009h . 

The bipolar wind of M82 galaxy has a co nical geometry 
with a total opening angle of ~ 30 degrees l|Ranalli et al.l 
2008). It has an X-ray luminosity of ~ 10 40 erg/s, temper- 
atures between ~ 10 and 10 s K (including the diffuse and 
warm material, dust and molecular gas) and a maximum ve- 
locity of the cold filaments (perpendicular to the disk) ~ 600 
km/s. 

The cold and hot gas metallicities in M82 are uncertain. 



Optical and infrared observations suggest for the cold gas a 
solar metal abundance, while the X-ray observations give 
an abundance for the hot gas ^ one-third the solar value 
^Strickland fc StevensllioOOr i. 

To understand the complexity of M82-like GWs many 
analytical and numerical two and three-dimensional studies 



fchevalier & Clegd 198d: Tenorio-Tagle & Munoz-Tunon 


19981: IStrickland & Stevens! l200d: 


Tenorio-Tagle et al. 


20031: [Rodriguez-Gonzalez et al.1 l008bl: 


CooDer et al.1 2008: 


Strickland & Heckmad 20091). IChevalier & Clegd (|l985l) 



were the first authors to propose a simple analytical model 
of a steady and adiabatic wind for M82. The same model 
was later revisited by several authors considering the 
radiati ve losses of the ejected gas with simpli fied assump- 
tions. iTenorio-Tagle fc Munoz-Tunod (| 19981 ) performed 
two-dimensional simulations of a biconical wind with 
radiative cooling, considering accreted matter from the 
external environment in order to determine the minimum 
mass of stars and gas i n the SB that was needed to offset 
the accreted material. Strickland fc Stevens! (|2000j) and, 
more recently, Strickland & Hec kmad (|2009h . helped by 
two-dimensional axi-symmetric simulations, studied the 
dependence of the wind dynamics, morphology and X-ray 
emission with the distribution of the ISM and the star 
formation history in the SB regions and studied also the 
distribution of the matter in dense clouds. Their results 
indicated that GWs, although efficient at transporting 
large amounts of energy to o utside of the galaxi e s, are 
inefficient at transporting mass. ITenorio-Tagle et al.l (|2003f) 
investigated also by means of two dimensional simulations 
the formation of filaments in a wind driven by SN e from 
several star clusters, while I Cooper et al l (|2008T l used 
three-dimensional (3D) simulations to study for the first 
time the formation of the wind in an inhomogeneous disk. 
They found that the presence of an inhomogeneous ISM in 
the disc may have an important effect on the morphology 
of the wind itself and on the distr i bution of the filaments. 
Finally, iRodriguez-Gonzalez et all (|008d lbl) also by means 
of 3D simulations have modeled the filaments in a GW 
as the result of the interaction between the winds from a 
distribution of SSCs and derived the condition necessary for 
producing a radiative interaction between the cluster winds, 
as in M82. In general, all these numerical studies indicate a 
coexistence in the wind of a low density, high temperature 
gas which expands freely and occupies a large volume with 
a cooler gas (with a filamentary structure which emits 
Ha lines) that expands more slowly and occupies a small 
fraction of the total volume. These models also predict a 
velocity for the hot phase (T ^ 10 8 K and n ^ 10 -4 cm -3 ) 
of ~ 1000 - 2240 km/s, and for the cold phase (T ^ 10 4 K 
and n < 10 2 cm -3 ) of ~ 600 km/s, in agreement with the 
observations. 

In this paper these previous analyses will be extended 
through a detailed numerical study of a SN-driven M82-like 
wind in the surrounds of the ejection region of the galaxy. 
We perform fully three-dimensional (3D) hydrodynamical 
simulations with radiative cooling in a gas in collisional 
ionisation equilibrium (CIE). The cooling function is com- 
puted taking into account the metal abundances both in 
the ISM and in the SNe ejecta. As in t h e previous studies of 
ITenorio-Tagle fc Munoz-Tunod l|l998h . ITenorio-Tagle et~afl 
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(2003) and iRodriguez-Gonzalez et all (|008bh . in our model 
the wind is generated by the ejected energy of the SNe 
hosted in the SSCs which are randomly distributed in time 
and space within a nuclear region of radius 150 pc. How- 
ever, in contrast to those studies, we employ a 3D galac- 
tic disk geometry with a multi-phase stratified and rotat- 
ing ISM which is particularly important to assess realis- 
tically the multi-phase complex wind's non-axisymmetric 
structures. Furthermore, with a more realistic treatment of 
the gas radiative cooling we will be able to track the forma- 
tion, evolution and stability of the dense and cold structures 
(filaments), which are formed as a result of the interactions 
between the dense radiatively cooled shells of the several 
supernova remnants (SNRs) generated in the SB, and coex- 
isting in equilibrium with the high temperature, low density 
wind component, as well as to determine more appropriately 
the fraction of the whole volume of the wind which is occu- 
pied by this cold phase. Also, the evolution of the SN ejected 
matter, (which is traced by a separated continuity equation 
as a passive scalar) will allow us to obtain for the first time 
in numerical simulations of GWs important information on 
the chemical evolution, as well as of the Ha and the soft 
X-ray emission of the wind, and the metal contamination of 
the intergalactic medium (IGM) due to the wind itself. 

In the following sections we will outline the main char- 
acteristics of our model (Sec. 2) and the results obtained 
from the numerical simulations, mostly highlighting the for- 
mation process of the filaments, the chemical evolution of 
the wind and the distribution of the most abundant species 
outside of the disk (Sec. 3). Finally, in Sec. 4 we will present 
a brief discussion and draw our conclusions, and in Sec. 5 a 
summary of the main results found in this work. 



our model we have $(r, z) = <E> star (r, z) + $>disk(r,z). We 
do not incorporate a dark matter halo component because 
we are interested in the behaviour of the wind evolving for 
only a short period of time (~ 10 Myr) and within a small 
volume in the surroundings of the galaxy, so that the grav- 
itatio nal effects of a dark matter halo are negligible in this 
case l|Strickland fc Steven j|200(f ). 

For an initially isothermal gas, the disk density distri- 
bution will therefore have the form: 



pd(r, z) = p d ,o exp 



$(r, z) - e 2 $(r, 0) - (1 - e 2 )$(0, 0) 



(4) 



where c„ 
and e 



quantifies the fraction of rotational sup- 



port of the ISM. To reduce the rotational velocity of the gas 
above the plane, we have assumed a simple model where the 
rotational support drops off exponentially with increasing 
height z and where the scale height fo r this reduction in the 
rotat ional velocity is z ro t = 5 kpc (see I Strickland fc Stevensl 
2000). Similarly, the distribution of the gas of the halo is 
assumed to be of the form: 



Ph(r, 



ph,o exp 



$(r, z) ~ e 2 $(r, 0) - (1 - e 2 )$(0, 0) 



(5) 



where c 2 T halo is the isothermal sound speed in the halo. 
We adopted a disk temperature, Tdisk ~ 6.7 x 10 4 K and 
a halo temperature Thaio ~ 6.7 x 10 6 K. Therefore, in each 
point of the system the total density is p(r,z) = pd{r, z) + 
Ph{r, z) and the total pressure is p(r, z) — pd(r, z)c 2 s T disk + 
ph{r, z)c\ T halo . The galaxy parameters adopted in this 
study are presented in Table [1] 



2 THE MODEL 

2.1 Setup for an M82-like starburst wind 

In our model, the gas of the disk is initially set in rotational 
equilibrium in the gravitational potential of the galaxy. The 
gravitational potential of the stars, <& a tar(r, z) is assumed to 
be generated by a stellar distribution with a spheroidal King 
profile 

P*,o 



p*(r,z) 



[1 + (r 2 + z 2 )/lu 2 ] 



2l3/2 ' 



(1) 



where p Ai o is the central density of the stars, and cjo = 
Ak /(ro + z 2 )=350 pc is the core radius. With this distribution 
the stellar gravitational potential is: 



<E>star(r, z) = 



GM S 



wo 



In 



j(vV 2 + z 2 /ujo) + y/l + (r 2 + I 
y/r 2 + z 2 /loo 

2 xlO 8 M© is the stellar mass 



2.2 Energy injection 

The SB activity in M82 is centred in the nuclear region 
with a diameter of ~ 500 pc. As described above, within 
this region there is a number of prominent and high surface 
bright ness clumps containing hundreds of y oung massive 
SSCs l|0'Connell et al.lll995l ; [Melo et alj|2005h . In our study 
we assume that the kp c-scale bipolar superwind observed 
in different wavelengths ( Shopbell & Bland-Hawthornl ll998l : 
IStevens fc Hartwelll [20031 : IStrickland fc Heckmanl 120071 ) is 
driven by the SN explosions occurring in the SSCs, and that 
the total luminosity injected into the ISM is ^^ ssc J Ji — 
10 42 erg s _1 , where Li is the luminosity of a single SSC 
and N S sc is the total number of SSCs. We also assume that 
Nssc ranges between a minimum of 10 and a maximum of 
K)Q and consider an inner nuclear region within a radius 



where M s tar=A'KpstarljJo- 

within the radius ujq. 

The gravitational potential of the gas disk is given by 
Miyamoto & Nagai model (1975), that is: 

GMdisk 



disk{r, z 



(3) 



+ ta + Vz 2 + b 2 ) 



where a = 222 pc and b = 75 pc are the radial and vertical 
scales of the disk, respectively, and Mdisk — 2x 10 9 Mq. In 



}/ w e(fi ,;b\ of 150 pc where this activity takes place. We further 
COTsraer that the total amount of SSCs may be formed in 
a sxar formation process along 1 or 10 Myr, depending on 
the model. Therefore, for each cluster we assume an aver- 
age luminosity Li = 10 42 /Af BS c erg s _1 and, since every SN 
injects E = 10 51 erg, we have JVsn,c = 10 42 (t&/ EoJVssc) = 
10 B /Nssc, where Msn,c is the total number of SNe per clus- 
ter, and where we assume that each single cluster is formed 
in an instantaneous burst which is characterized by SN ac- 
tivity occurring over a time t& = 30 Myr (i.e. the lifetime of 
the less massive SN II prog enitor star) and by a SN heat- 
ing e fficiency HE = 100% (|Melioli fc de Gouveia Dal Pinol 
120041 ). where the heating efficiency, HE, corresponds to the 



4 C. Melioli, E. M. de Gouveia Dal Pino, F. G. Geraissate, 

Table 1. Parameters for the galaxy setup obtained from M82 
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Figure 1. Number of SN explosions in the first 10 Myr for each 
model. Solid (black) line: SC40-10; dashed (blue) line: SC10-01; 
long-dashed (green) line: SC40-01; dotted (red) line: SC10-10 and 
dotted-dashed (magenta) line: SC100-10 



fraction of the SN explosion energy that remains effec- 
tively stored in the ISM gas (as kinetic and thermal en- 
ergies) and is not radiated awayQ Therefore, an SSC in- 
jects mass and energy at rates M = M e j W 42 / (EoAfssc) and 
L w = 10 42 /A 'ssc, respectively, where M e j= 16 Mq is the 
mean mass rele a sed b y a single SN explosion. According to 
iMarcolini et al.l (|2007l ). we also assume that each SN explo- 
sion injects 3 Mq of metals into the ISM. Finally, to set 
the clusters spatially and temporally within the core of the 
galaxy, we associate randomly to each i-th cluster a position 
P l = (n,<f)i), where < n < 150 pc and < < 2ir, and 
a time t* in the range < f < 1 Myr or < t* < 10 Myr 
depending on the model. 

A summary of all the parameters adopted for the inter- 
mittent wind models studied here are given in Table [2] and 
Figure [1] shows histograms representing the number of SN 
explosions during the first 10 Myr for each model. 



2.3 Numerical methodology 

To simulate the evolution of the superwind driven by SN ex- 
plosions in an M82-like SB galaxy we use a modified version 
of the numerical adaptive mesh refinement hydrodynamica l 
code YGUAZU (|Raga et alj|200d . l2002t rMelioli et al.ll2008l . 
2009) that integrates the gas-dynamic equations with the 
flux vector splitting algorithm of Van Leer (1982). It also 
includes a parametrized cooling function in the energy equa- 
tion that allows the gas to cool from ~ 10 7 ' 5 to 10 4 K with 
errors smaller than 10% and which is calculated implicitly 
in each time step for each grid position. The 3D binary, hi- 
erarchical computational grid is structured with a base grid 
and with a number of nested grids whose resolution doubles 
going from one level to the next one. We have run the mod- 
els with a maximum resolution of 5.8 pc per cell. For one 
model (SC40-01), which we take as the reference model, we 
have adopted five grid levels covering a box with physical 
dimensions of 1.5x1.5x3.0 kpc in the x, y and z directions, 
respectively, while for all the remaining models we have con- 
sidered a box with physical dimensions of 1.5x1.5x1.5 kpc. 
The maximum dimensions of the computational domain for 
each model are also presented in Table 2. The energy associ- 
ated to each SSC is injected as thermal energy in a single cell 
of the highest grid level of the box. In order to account for 
the metal abundance injected by the SN explosions into the 
system, we have introduced an extra separated continuity 
equation into the code in the form dp z /dt + V • (p z v) = S z , 
where p z is the metal density, and S z is the SN metal source. 
This also allows us to easily track the metallicity evolution 
in the wind. 

Finally, for computing the radiative cooling we em- 
ployed a cooling function A(T) considering an optically thin 
gas in the ionization equilibrium regime, which implies 

C(p,T) = p 2 A{T)Z erg cm" 3 s _1 (6) 

llMcWhirter et al.1 Il975l; lAldrovandi fc PequignotJ Il973l ; 



Ide Gouveia dal Pino fc Bend ll993af T where Z is the total 
gas abundance computed at each time and at each point of 
the grid. 



3 RESULTS 

In this section we present the results for the different models 
listed in Table [2] exploring a set of SSC parameters aiming 
mostly at the examination of the injected energy, the multi- 
phase ISM and the chemical evolution of the system. All the 



We note that due to the short lifetime of the starburst, almost 
all the injected energy in the wind is due to SNe II only. For this 
reason, in this study we neglect the contributions due to type I 
SNe explosions. 



2 We have also performed several tests with a larger resolution 
of 2.9 pc per sell for 2 Myr and found that the results were very 
similar to those with a 5.8 pc resolution per cell. For this reason, 
we adopted the latter scheme for the study presented here as it 
allowed us to save computation time without loosing information 
and quality in the results. 
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Table 2. Wind parameters used in our hydrodynamic simulations for 5 different models. For each Model: 
a) Number od SSCs; b) Duration of the SF process (in Myr); c) Injected luminosity per SSC (in erg/s); d) 
Injected matter per SSC (in MQ/yr); e) Total energy injected after the first 5 Myr (in units of 10 55 erg); 
and /) Physical dimensions of the box (in kpc). 



Model SSCs (J) SF (Myr) I7/SSC (erg/s) Af/SSC (M /yr) E 5 (10 55 erg) Box Size (kpc) 
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models, except the reference one (SC40-01), were evolved for 
5 Myr right after the first stars in an SSC (with masses 120 
M©) have begun to explode, which is an appropriate time to 
investigate the global properties of the wind build up, from 
the first SN explosions up to a nearly steady state phase 
reached after ~ 3.5 Myr. The reference model, instead, was 
evolved for 8.5 Myr, in order to follow the gas outflow evolu- 
tion up to a height of 3 kpc. The results here obtained may 
be directly compared with observations of GWs in order to 
determine, for instance, the best set of initial conditions able 
to drive their outflows. 

Table [3] summarizes the global energy budget for all the 
models after the reference times t=3.5 and t=5 Myr. These 
energies will help us to understand the wind evolution in the 
discussion below. 



3.1 Model SC10-10 

The first model (SC10-10) considers 10 SSCs formed along 
10 Myr, each of them having a total stellar mass of 10 7 Mq, 
a total SNe number of 10 5 and a mean separation of ~ 50 
pc. With such configuration, the SN explosions occurring in 
the first SSCs have more than sufficient energy to drive an 
outward propagating shock wave which sweeps up the ISM 
very quickly. In fact, as a single SN explosion generates a 
supernova remnant (SNR), the whole set of SNe explosions 
due to the first SSCs generates several SNRs which interact 
with each other still in the nuclear region of the disk and 
drive an outward propagating large scale shock front. This 
sweeps up the ISM into a thin and dense shell that envelopes 
a large, hot, low-density cavity. Such structure, usually de- 
nominated superbubble, has been also detect ed in former nu- 
merical studies of GW formatio n (see, e.g.. iTomisaka et al.l 
Il98ll : iTomisaka fc Ikeuchill 19861 ). It is maintained at a high 
temperature and a low density due to the continuous in- 
jection of SN energy. New superbubbles generated by later 
SSCs will expand into such a rarefied ISM without form- 
ing new dense shells and thus suffering negligible radiative 
losses. Therefore, with the setup of model SC10-10, almost 
all the SN energy is transferred to the ISM without being 
radiated away and the result is a highly energized super- 
bubble driven by the energy ejected from each SSC, with- 
out forming a wind with filamentary structures immersed in 
it. Nonetheless, at least the shock compressed external shell 
that develops in the beginning is radiativelly cooled and suf- 
fers fragmentation due to the combination of the develop- 
ment of the Kelvin-Helmholtz (K-H) and Rayleigh- Taylor 
(R-T) instabilities as the dense gas of the shell is pushed 



into the rarefied gas of the disk halo (e.g. lMelioli et al1l2005l ; 
Ide Gouveia dal Pino fc Benzlll993bl ). 

This is the resulting scenario that we see in Figure [2] af- 
ter 3.8 Myr when the upper part of the wind has already left 
the computational domain. The top-left panel shows the gas 
column density and we note that the only dense structures 
are formed along the expanding shocked shell. No filaments 
or clouds are formed within the wind flow and consequently 
the gas of the disk carried out by the wind and lost to the 
IGM is essentially the interstellar gas swept by the first gen- 
eration of SN explosions. This is also suggested by the top- 
right and bottom- left diagrams of Figure [2] The low density 
phase of the wind, which is enriched by the metals ejected 
by the SNe, has a high abundance ~ 10 times larger than 
the solar abundance. Nonetheless, most of the metals remain 
in the nuclear region of the galaxy. In this case we also note 
that the thermal energy radiated by the wind comes from 
the external shell, characterized by a temperature of ~ 10 4 
K (bottom-right panel of Figure [5J). 

Therefore, the results above suggest that when the gas 
outflow is generated by only a few very massive SSCs most 
of the soft and/or hard X-ray emission will be detected only 
behind the shocked regions of the shell-like structure around 
the GW, provided that the reverse shock where the wind 
impacts the outward propagating shock is right behind the 
later. 

Investigating carefully the features of the gas, we real- 
ize that in model SC10-10 there are two well separated gas 
phases. The first one associated to the internal part of the 
superbubble which has a density smaller than 10 _1 cm -3 , 
a temperature of ~ 10 7 K and a velocity between 100 and 
2000 km s~ . The other gas phase is associated to the exter- 
nal shell of the superbubble, with a density larger than 10 _1 
cm~ 3 , a temperature of ~ 10 4 K and a velocity normal to 
the disk direction of ~ 1500 km s . These features can be 
also distinguished in Figure|3] where we have plotted the gas 
mass as function of its vertical velocity, for three different 
ranges of gas densities, and in Figure[4] where the total mass 
of each phase with a vertical velocity ^ 50 km s _1 is plot- 
ted as function of z (the height above the disk) . We observe 
that the intermediate-phase, characterized by a density of 
~ io _1 -io- 2 cm , has a similar mass to the high-density 
phase, of ~ 3xl0 5 Mq. Its velocity follows the velocity dis- 
tribution of the high-density phase up to ~ 500 km s -1 and 
then, follows the distribution of the low density phase. This 
means that its evolution does not occur separately from the 
other phases. In the disk (where the velocities are smaller), 
this phase is related to the high density phase, while in the 
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Table 3. Energy budget (in units of 10 55 erg) at t=3.5 Myr (-E3.5) and t=5 Myr (E 5 ). For each Model: 
a) energy injected by the SN explosions; b) total energy stored into the system; c) total kinetic energy; d) 
vertical (or poloidal) kinetic energy; e) total thermal energy; /) total energy lost by radiative cooling; g) 
fraction of the injected SN energy which is radiated away. 



Model 


Esjv,3.5 


Esj/,3.5 


E fc tot ,3.5 


E fcj,3.5 


Eth,3.5 


Eiost,3.5 


E r ad,3.5 (%) 


SClO-10 


2.20 


1.10 


1.28 


0.45 


0.20 


1.50 


68 


SClO-01 


5.60 


2.35 


2.00 


0.98 


0.33 


4.10 


73 


SC40-10 


0.86 


1.18 


0.98 


0.20 


0.19 


0.53 


62 


SC40-01 


4.90 


2.80 


2.46 


1.50 


0.39 


2.95 


60 


SClOO-10 


1.10 


1.14 


0.95 


0.10 


0.23 


0.81 


74 


Model 


EsJV,5 


Esj/,5 


E fc tot ,5 






E(ost,5 


E ra d,5 (%) 


SC10-10 


3.20 


1.75 


1.49 


0.42 


0.26 


2.30 


72 


SC10-01 


7.30 


2.15 


1.90 


1.09 


0.30 


6.00 


82 


SC40-10 


1.60 


1.25 


1.08 


0.26 


0.17 


1.20 


75 


SC40-01 


7.20 


2.70 


2.35 


1.30 


0.35 


5.35 


74 


SC100-10 


2.10 


1.30 


1.04 


0.28 


0.26 


1.65 


79 



wind (where the velocities are larger), it is connected to the 
low density phase. Figure [3] also highlights the presence of a 
new generation of accelerated shells which is represented by 
the peak velocities between 1200 and 1700 km s _1 . It is well 
separated from the high density gas, which has velocities be- 
tween 50 and 900 km s" 1 and is probably distributed along 
the external shell pushed by the outflow. The total mass of 
this new generation of shells is small, because almost all the 
gas of the nuclear region has been already swept by the first 
generation of SN explosions and therefore, no filaments or 
structures are driven along the wind. 

The results above indicate that the SN-driven gas out- 
flow from the nuclear region of a galaxy with only a few 
highly massive SSCs radiates about 68% of the total energy 
injected by the SNe. In fact, after 3.5 Myr, when the wind 
has already filled completely the physical domain of our sim- 
ulation, the sum of the kinetic and thermal energies of the 
system is ~ 1.5 xlO 55 erg, i.e. almost twice the energy of 
the system at t = 0, corresponding to 8.5 x 10 54 erg, while 
the total energy injected by the SN explosions is about 2.2 
x 10 55 erg. Therefore, nearly 2/3 of the SN energy is radiated 
away (mainly in the nuclear region and in the external shell), 
and the remaining 1/3 goes to the kinetic energy (20% of 
the total injected energy) and thermal energy (10% of the 
total injected energy) of the gas. This result is consistent 
with the evolution of the multi-phase ISM shown in Figure 
[4] Over the simulated period, the amount of gas with a cool- 
ing time shorter than the dynamical time, represents ~ 60% 
of the total gas flow and thus assuming that the SN energy 
is spread approximately uniformly, it is reasonable to be- 
lieve that all the energy in the shells is radiated away, while 
the remaining energy in the low and intermediate-density 
gas phases increases the gas pressure and consequently, the 
thermal and kinetic energy of the GW. 



Column Density 



-0.4 0.0 m 
Weighted Abundance 



L ] 




Figure 2. Model SC10-10: color-scale map of the gas column 
density (top-left panel), the metal column density (top-right), the 
total weighted abundance, (i.e. ^y™g X Pz{%> *)/ X/y— o* ^^ x ' 
where pz(x,z) and p(x,z) are the metal and total gas density 
in the space coordinates x,y, respectively) ( bottom-left) and the 
total emission (bottom-right) at at a time t=3.8 Myr. The column 
density is shown in units of cm -2 , the abundance is normalized 
by the solar abundance and the emissivity is in erg s" 1 . All the 
quantities are expressed in log-scale. In this model, the SSCs have 
a mass of 10 7 Mq, a SN number of 10 5 and are formed in a 
starburst process of 10 Myr. 



3.2 Model SC10-01 

The second Model (SC10-01) considers, as in the previous 
one, 10 SSCs with a stellar mass of 10 7 Mq, a SNe number 
of 10 5 and a mean separation of ~ 50 pc. However, in this 
case we assume that all the super clusters are formed along 
1 Myr only and therefore, in contrast with the first model, 



every SSC develops a hot bubble surrounded by a thin dense 
shell at about the same epoch. The SN explosions occurring 
in the SSCs have sufficient energy to drive the bubble ex- 
pansion at the same time that the collisions between the 
shells lead to the formation of an interconnected network of 
dense filaments permeated by tunnels containing very hot, 
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Figure 3. Model SC10-10: Gas mass distribution as function of 
the vertical velocity for three different phases of the gas at t=3.8 
Myr. Solid line: gas with number density larger than 1 cm -3 ; 
dashed line: gas with number density between 10~ 2 and 1 cm -3 ; 
dotted-dashed line: gas with number density smaller than 10" 2 
cm -3 . The velocity is expressed in km s" 1 and the mass is in 
units of solar mass (in log-scale). In this model the SSCs were 
formed over a period of 10 Myr. 



Figure 4. Model SC10-10: Total gas mass distribution as function 
of z (height above the disk) at three different times: t=1.5 Myr 
(dashed line), t=2.5 Myr (dotted-dashed line) and t=3.8 Myr 
(solid line). Top- left panel: gas with number density larger than 
1 cm -3 ; top-right panel: gas with number density between 10" 2 
and 1 cm -3 ; bottom-left panel: gas with number density smaller 
than 10 — 2 cm -3 . Bottom-right panel: comparison between the 
total gas mass at t=3.8 Myr in the high (solid line), intermediate 
(dashed line) and low (dotted-dashed line) density phases of the 
gas. Mass is expressed in Mq (log-scale) and z is in kpc. 



tenuous ISM. In fact the SNRs generated by the SN explo- 
sions propagate into an ambient medium with a density > 
1 cm -3 , and e nter in the rad iative phase in a time shorter 



than 3x l0 4 vr jMcCravl 1987). that is, at radii smaller than 



-20 pc (|Melioli et alj|2006t ). Therefore, their external shells 
are expected to become very dense (n a h/n ^ 10) and cold 
(T ^ 10 4 K) before interacting with each other. Collisions 
among them lead to the formation of very dense filaments 
(up to 1000 times denser than the hot gas) and dense sur- 
faces (up to 100 times denser than the hot gas) in times 
between 10 5 and 10 6 yr according to our simulation. R-T 
and K-H instabilities are also able to fragment these struc- 
tures and produce new filaments, in a process that can con- 
tinue during all the SN activity, until th e average gas den- 
sity of the SB is ^ 10" 1 cm" 3 (see e.g.. lMcCravlll987l and 
iMelioli et al.|[200(3 ). In spite of the initial formation of fila- 
ments and dense structures, the superbubble resulting from 
the interaction of the individual SSC bubbles evolves as in 
Model SC10-10 and when also the external shell fragments 
due to R-T and K-H instabilities, the hot gas is able to es- 
cape from the galaxy without carrying very large amounts of 
metals or disk gas. Nonetheless, Figure [5] which maps the 
column density edge-on distributions of each gas phase in 
different velocity ranges, clearly shows that a bunch of fila- 
ments and a large amount of metals launched together with 
the hot and low density gas (top-left and top-right panels 
of Figure [5| and consequently, a large fraction of energy is 
radiated away by the filaments themselves. 



3.3 Model SC40-10 

Now we are going to analyse the results of our simulations 
considering a SB region with intermittent injection of mass 
and energy from a larger number of SSCs (40) having a 
stellar mass of 2.5 x 10 6 M , a SNe number of 2.5 x 10 4 and 
a mean separation of ~ 40 pc. In this case we assume, as in 
Model SC10-10, that the SSCs are formed over a time of 10 
Myr (Model SC40-10). 

The resulting gas outflow is very similar to the one ob- 
tained in Model SC10-10 but it has an amount of filaments 
up to a height ~ 600 pc. Also in this case, despite the high 
metal abundance of the wind flow (~ 10 Zq), most of the 
metals ejected by the SNe remain in the disk, and only a 
small fraction is transported to outside of the system to- 
gether with the large external shell formed by the activity of 
the first SSCs. Finally, the gas emission of the central region 
(up to 500 pc) is about 1000 times larger than the emission 
of the wind and 100 times larger than the emission of the 
external shell. Therefore, in such a case, we expect an X-ray 
emission coming only from the SB core of the galaxy and a 
weak radio and infra-red emission from the shell fragments 
at ~ 1.5 kpc from the disk, characterized by a maximum 
temperature of ~ 10 4 K. 

Although the main features of the wind are different 
from those observed in M82, we note that an increase in 
the number of SSCs from model SC10-10 to model SC40-10 
favours the formation of a multi-phase gas in the latter. The 
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Figure 5. Model SC10-01: map of the gas column density for 
different density and velocity ranges at a time of 3 Myr. top- 
left panel: n ^ 0.5 cm -3 and v ^ 50 km s — 1 ; top-right panel: 
0.01 n ^ 0.5 cm -3 and v ^ 50 km s^ 1 ; bottom-the left panel: 
n ^ 0.01 cm -3 and v ^ 50 km s — 1 ; bottom-right panel: n ^ 0.5 
cm -3 and v ^ 300 km s _1 . The column density is expressed in 
cm -2 . The "white" portions in each panel represent regions with 
total column density smaller than 8xl0 17 cm -2 . 

dense, cold gas phase shows a peak between velocities of 600 
and 800 km s _1 , corresponding to a total gas mass of about 
3000 Mq in accelerated filaments in the nuclear region. At 
the same time, both the intermediate and low density gas 
are accelerated up to velocities of 2000 and 2300 km s -1 , re- 
spectively, and the intermediate (warm)-phase velocity dis- 
tribution overlaps with that of the denser (cold)-phase until 
800 km s _1 , which corresponds to the maximum velocity of 
the filaments. 

To understand better the multi-phase gas distribution, 
as before, we have analyzed the edge-on column density of 
each gas phase, which shows that in the first 500 pc the cold 
(high density), warm (intermediate density), and hot (low 
density) gas phases are completely mixed up. Above 500 pc, 
instead, the phases are decoupled, and no clouds or filaments 
are observed together with the low density launched flow. 

In this model the energy radiated away is about 75% of 
the total energy injected by the SNs. After 5 Myr the sum 
of the kinetic and thermal energies of the system is ~ 1.25 
xlO erg, while the total injected energy is ~ 1.6 xlO 
erg. This means that about 3/4 of the SN energy is lost by 
radiative cooling, while the remaining 1/4 goes into the ki- 
netic energy (15% of the total injected energy) and thermal 
energy (10% of the total injected energy) of the gas. Also in 
this case, the energy distribution is consistent with the com- 
position of the multi-phase ISM. After 4 Myr, the amount 
of gas in the intermediate and high density phase represents 
~ 65% of the total gas flow, but we note that earlier this 
fraction was larger than 80%. This means that the energy 



was radiated away with high efficiency and proportionally to 
the total amount of gas with density larger than 0.01 cm -3 . 

3.4 SC40-01 - the reference Model 

The fourth Model (SC40-01) considers 40 SSCs (as in Model 
SC40-10) with a stellar mass of 2.5 x 10 6 M Q , a SNe number 
of 2.5 x 10 4 , a mean separation of ~ 40 pc, but with the for- 
mation of the SCCs over 1 Myr only. This model, which we 
denominate our Reference Model, due to the larger energy 
power injected (see Table 3) is able to produce and launch 
a multi-phase wind to greater distances. For this reason, we 
have run this simulation in a larger box domain with phys- 
ical dimensions of 1.5 x 1.5 x 3 kpc. 

As described in Model SC-10-01, the larger number of 
SSCs and the shorter time (1 Myr) at which they begin 
to be active (which lead to a larger injected energy power), 
have created the conditions to produce an environment with 
rich filamentary structure and with a thermal pressure large 
enough to drive both the low and high density phases of 
the wind very quickly to large distances. In a few Myr, a 
multi-phase ISM is launched in the IGM and a GW with a 
large amount of dense structures embedded in a hot low den- 
sity gas develops. Figure [6] shows this situation. Above the 
disk, at z 3> 300 pc, the edge-on column density distribution 
shows the presence of a network of elongated filaments in the 
direction of the rarefied flow. Due to their high density (in 
spite of their low abundance), the filaments carry with them 
a large amount of metals and therefore, in this case, repre- 
sent an efficient mechanism to transport the chemical species 
ejected by the SN explosions above the disk. Also, the metal 
abundance in the wind is large, between 5 and 10 Z@, and 
these values are in agree ment with those obtained in previ- 
ous studies of M82 wind (|Strickland fc Heckmail 2009). The 
launched multi-phase medium radiates energy at all heights 
(see Figure[Bl bottom-right panel) and thus may be observed 
both in X and radio wavelengths. 

If we look at the velocity distribution of the different 
gas phases (Figure P7J), we note that gas at high density 
with velocities between 150 and 800 km s _1 has mass al- 
most 10 times larger than in previous models. This means 
that the dense (cold)-phase is being launched from the SB 
region along with the hot phase in the wind. This trend 
is confirmed by the multi-phase evolution shown in Figure 
[5] The dense, intermediate and low density gas phases are 
completely superimposed, each one with its own velocity. 
In particular, the high velocity clumps and fragments are 
spread and embedded all over the wind. The total mass of 
the high density phase in the flow is ~ 4 x 10 6 Mq, while 
the gas of intermediate density has a mass of ~ 1.5 x 10 6 
Mq, and the gas of low density has ~ 5 x 10 Mq (Figure 
[9|. The high density gas in the wind is therefore, ~ 2/3 of 
the total GW mass, while only 8% of the GW mass is due 
to the hot and high-velocity, low density gas. Due to the 
large fraction of high density gas in the flow, most of the 
injected energy is lost by radiative cooling. After 5 Myr the 
total energy of the system is 2.7 xlO 55 erg, while the total 
injected SN energy is 7.2 xlO 55 erg. Since the initial energy 
is 8.8 xlO 54 erg, this implies that the energy radiated away 
by radiative cooling is ~ 5.35 x 10 55 erg, that is, ~ 74% of 
the injected energy. 

As the main features of the GW in M82 seem to be sim- 
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Figure 6. Model SC40-01: the same as in Fig. [2] at t=3.5 Myr, 
except that here we consider 40 SSCs forming over 1 Myr with a 
mass of 2.5 xlO 6 M Q and a SN number of 2.5 xlO 4 . 



Figure 7. Model SC40-01: the same as in Fig. [3] except that here 
we consider 40 SSCs forming over 1 Myr with a mass of 2.5 XlO 6 
M Q and a SN number of 2.5 XlO 4 . 



wind. These values are similar to those inferred fr om recent 
observations (see, e.g., IWestmoquette et al. I l009al lbl). More- 
over, looking at the vertical velocity distribution □ we see 
that the filaments produced in the core move in the wind 
with a velocity between ~ 100 and 800 km s _1 , while the 
gas with high thermal pressure that drives the superwind 
moves in different streams and channels with a typical veloc- 
ity of ~ 2000 km s" 1 . Considering th at the escape velocity of 
the M 82 wind is ~ 400 km s" 1 (e.g.. lStrickland fc Heckmanl 
2009), we can conclude that only the low density gas and 
a little fraction of the filaments may really escape from the 
galaxy. 

Finally, in Figure [11] we show the average abundances 
for each gas-phase of the flow. The faster and smaller density 
phase, with a velocity between 500 and 2000 km s _1 , is 
characterized by a high metal abundance, ~ ten times larger 
than the solar abundance. On the other hand, the slower 
and denser phase, with a velocity between 100 and 800 km 
s" 1 , representing ~ 70% of the total flow mass, has a metal 
abundance between 1 and 4 times the solar value. 



ilar to those obtained in this Reference Model (Sc40-01), we 
can investigate in more detail the physical conditions of the 
energized gas in this model. Figure [TO] depicts the face-on 
gas density, temperature, pressure and vertical-velocity dis- 
tributions at an earlier time t=2.5 Myr, at a height of 50 pc 
above the disk where the coexistence of different phases is 
better highlighted. A straightforward feature in this figure 
is the presence of multi-phase structures with high temper- 
ature (T ^ 10 7 K), pressure (p/k ~ 4 x 10 7 ) and density 
(n~5x 10 2 cm -3 ) coexisting with an underlying high ve- 
locity (v ~ 2000 km s" 1 ), low density (n sC 10~ 2 cm" 3 ) 



3.5 Model SC100-10 

The last Model (SC100-10) considers 100 SSCs built up over 
10 Myr with a stellar mass of 10 6 M©, a SNe number of 
10 4 , and a mean separation of ~ 25 pc. In this case, as in 
the models SC10-10, SC10-01 and SC-40-10, the numerical 
simulation was performed in a box domain with dimensions 
1.5 x 1.5 x 1.5 kpc. 

The results are similar to those of model SC40-01. The 



3 We have examined several cuts along and in the normal direc- 
tions of the flow which are not presented here. 
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Figure 8. Model SC40-01: map of the gas column density for 
different density and velocity ranges, at a time of 4.5 Myr. Top- 
left panel: n ^ 0.5 cm -3 and v ^ 50 km s — 1 ; top-right panel: 
O.Of ^ n ^ 0.5 cm -3 and v ^ 50 km s — 1 ; bottom-left panel: 
n H 0.01 cm -3 and v ^ 50 km s — 1 ; bottom-right panel: n ^ 0.5 

cm -3 and v ^ 300 km s _1 . The column density is expressed in 
™-2 



larger number of SSCs combined with a smaller power in- 
jected per SSC (which are formed for longer time) generates 
an ambient which is rich in filaments and with a high ther- 
mal pressure that accelerates the filaments and clouds to 
high velocities together with the hot wind flow. Also in this 
case the multi-phase ambient with dense structures embed- 
ded in the hot low density gas provides alone 70% of the 
total flow mass (Figure [12)) . 



Figure 9. Model SC40-01: the same as in Fig. [4] except that here 
we consider 40 SSCs built up over 1 Myr with a mass of 2.5 xlO 6 
Mq and a SN number of 2.5 XlO 4 at a time of 1.5 Myr (dashed 
line), 3.3 Myr (dotted-dashed line) and t=4 Myr (solid line). 
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Figure 10. Model SC40-01: Color-scale map of the face-on M82 
gas density (top- left panel), temperature (top-right), vertical- 
velocity (bottom-left) and pressure (bottom-right) distributions, 
at a height of 50 pc above the disk and at a time t=2.5 Myr. The 
density is shown in units of g cm -3 , the temperature in K, the 
pressure in dyne cm -2 and the velocity in km s — 1 . 
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Figure 11. Model SC40-01: map of the average metal abundance 
for different density and velocity ranges, at a time of 4.5 Myr. 
Top-left panel: n > 0.5 cm -3 and v ^ 50 km s^ 1 ; top-right 
panel: 0.01 ^ n ^ 0.5 cm -3 and v ^ 50 km s -1 ; bottom-left 
panel: n ^ 0.01 cm -3 and v ^ 50 km s — 1 ; bottom-right panel: 
n 0.5 cm - 3 and v ^ 300 km s — 1 . The abundance is expressed 
in unit of solar abundance. 



The main difference with regard to the previous model 
(SC40-01) is the process of clump formation. In model SC40- 
01, after 2 Myr the filamentary structure has a total mass of 
~ 6 x 10 6 M , while after 5 Myr it decreases to ~ 3.5 x 10 6 
M©. In model SC100-10, on the other hand, the filamentary 
mass increases with time from 6 x 10 s Mq at t = 2 Myr 
to ~ 2.3 x 10 6 M Q at 5 Myr. Around 4 Myr the total mass 
of the filaments is about the same in both models, but in 



Figure 12. Model SC100-10: the same as in Fig. [I] except that 
here we consider 100 SSCs with a mass of 10 6 M Q and a SN 
number of 10 4 at a time of 2.5 Myr (dashed line), 3.5 Myr (dotted- 
dashed line) and t=5 Myr (solid line). 



model SC40-01, a large amount is quickly produced and then 
partially ablated by the wind flow, while in the second case 
(model SC100-10), the filaments are produced more slowly 
as the superbubbles expand and collide. This different evo- 
lution of the denser gas phase affects the radiative cooling 
too. For model SC100-10, after 5 Myr the total energy of 
the system is 1.3 xlO erg, while the injected SN energy 
is 2.1 xlO 55 . This means that the energy radiated away is 
~ 1.65 x 10 55 erg, or ~ 80% of the total injected energy, 
therefore, a little larger than that in model SC40-01. 



4 DISCUSSION AND CONCLUSIONS 

In this study, we have investigated the formation of galac- 
tic winds (GWs) driven by SNe explosions, with particular 
focus on the M82 galaxy. 

Several driving processes of GWs have been pro- 
posed over the last years, including the r mal pressure 
due to supernova heating jLarson 19741; foekel fc Silk 
19861 ; iMac Low fc Ferraral Il999l; b'Ercole fc Brighenti 



19991; IStrickland fc Stevens! |2000|; iMelioli et al.ll200Sl. 120091; 



Hill et al.ll2012ri, radi a tion p r essure on dus t grain s (|Martinl 



20051 ; iMurrav etalJ 120051; iNath fc Silkl I2009T) 



ray p ressure ( Breitschwerdt et al.l Il99ll : lEverett et al.l 
I2OOSI; lUhTig et all l2012t) an d supersonic turbulence 
1 Scannapieco fc Briiggenl |2010| ). In most of these pro- 
cesses, SNe are the main source of energy injection which 
is then reprocessed. In the particular case of SB galaxies, 
where the SFR is high, we expect that the SN energy 
dominates over all the other processes and therefore, in 
the present study we have only considered the role that 
SN shock fronts play on the evolution of the ISM (at 
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small scales) and on the GW evolution (at large scales). 
In the central SB of M82 there is a large number of super 
stellar clusters (SSCs) which are permeated by a complex 
distribution of gas with different densities, pressures and 
temperatures. As expected, most of our 3D radiative 
cooling hydrodynamical numerical simulations generated 
outflows with a complex morphology which captured 
the main features of an M82-like wind, i.e., a reach cold 
filamentary structure embedded in a hot low dense gas (see 
model SC-10-01, SC-40-01 and SC-100-10), starting right 
in the nuclear region and therefore in agreemen t with the 
observations (see, e.g. IWestmoquette et al.ll009al ). 

Although previous works have already studied the M82 
wind taking into account the effects of radiative losses, here 
we have performed a detailed analysis exploring the effects 
of different initial conditions in order to assess the connec- 
tion between the characteristics of the SB environment and 
the main features of an M82-like wind. More specifically, we 
tested the effects of two relevant parameters in the evolution 
of a SB GW: the total number of SSCs which are randomly 
distributed in space and time in the SB nuclear region (and 
their stellar mass) and the rate of star formation (or in other 
words, the build up rate of the SSCs), which ultimately trace 
the evolution of the SNe injected power in the wind (see be- 
low). 

In the following we discuss and summarize the main 
results obtained, particularly emphasizing the evolution and 
characteristics of the multi-phase environment, the energy 
budget of the system, the GW metal abundance, and the 
mass and metals lost by the galaxy to the IGM. 

4.1 Multi-phase wind environment 

To study the formation and evolution of a multi-phase envi- 
ronment we have considered three typical gas phases which 
are mainly distinguished by their density. We denominated 
high-density phase the gas with a number density n ^ 0.5 
cm -3 ; intermediate-density phase the gas with a density be- 
tween 10 -2 and 0.5 cm -3 ; and low-density phase the gas 
with a density n 10~ 2 cm' 3 . Considering in all models 
the same final total injected SN luminosity of 10 42 erg/s, 
but with different injection rates over the first 10 Myr (see 
Figure [TJ, we have shown that GWs are build up with a 
low-density gas component with typical velocities between 
50 and 2200 km s _1 and with a high-density gas component 
with velocities between 50 and 800 km s _1 . 

The SB region contains a total amount of gas of ~ 
3xl0 7 Mo, but only a fraction of this gas is accelerated 
by the SN explosions and driven into the wind flow, while 
the remaining part is pushed to the outskirts of the SB en- 
vironment (i.e. to the external regions of the disk) and part 
remains among the SSCs in a turbulent state. 

Because in all models the SB region is very quickly filled 
with the superbubbles produced by the first SSCs, the SB 
mass injected in the wind is nearly the same in all mod- 
els, ~ 65% of the total mass at the beginning of the SN 
activity, i.e. ~ 2 Myr, settling then to about 75% at the 
end of the simulations, at 5 Myr. However, depending on 
the model, the diffuse gas may be either converted in dense 
filaments and clumps due to the compression and radia- 
tive cooling of the swept material by the SN shock waves, 
or just ejected from the SB region in a hot, low density 



wind surrounded by a thi n, dense shell. This i s cons istent 
with the results found by iTenorio-Tagle et al.l (120031 ) and 
iRodrfguez-Gonzalez et all l|008al lbf). In both cases a GW de- 
velops, but with distinct morphology. In the first case we 
have "filamentary- like" GWs, as observed in M82 and also 
in models SC10-01, SC40-01 and SC100-10, while in the sec- 
ond case we have "supperbubble-like" GWs where the su- 
perbubble expands until its external super-shell fragments 
allowing the hot gas to flow into the IGM and where the 
only structures are those generated by the fragmentation of 
the super-shell, as in models SC10-10 and SC40-10. We note 
that, regardless of the total power injected which is the same 
for all models (10 42 erg/s), the spatial and temporal distri- 
bution of the SSCs within the SB region and thus the first 
million years of SNe activity (Figure[T}, are determinant fac- 
tors to allow the formation of a rich filamentary structure 
in the wind. If the SSCs form fast enough and if they are 
in a sufficiently large number then, several SN shells evolve 
overall the region and can interact nearly simultaneously, 
building up a rich filamentary structure which is dragged by 
the wind. 

Considering the discussion above, we may try to estab- 
lish the physical conditions which are necessary to form a 
multi-phase ISM and an M82-like GW. We have found that 
a larger amount of clumps and filaments is mainly generated 
during the first interactions of the SNe superbubbles when 
their volume filling factor is // ^ 1 and a large fraction 
of unperturbed gas is still in the nuclear region. Therefore, 
it is the number of active SSCs before // becomes //=1 
that determines whether the wind will be a "filamentary" 
or a "superbubble-like" wind. Our results indicate that // 
becomes equal to unity at £=1 Myr for model SC10-10, at 
£=0.6 Myr for model SC10-01, at £=1.3 Myr for model SC- 
40-10, at £=0.65 Myr for model SC40-01, and at £=1.9 Myr 
for model SC-100-10. At these times, the number of active 
SSCs is 3 for SC10-10, 7 for SC10-01, 4 for SC40-10, 26 for 
SC40-01 and 23 for SC100-10. This means that the larger 
the number of SSCs formed before // becomes // = 1, the 
more complex and rich the wind filamentary structure will 
be. When the number of active SSCs is ^ 7-10 and // ^1 
(which corresponds to a mean separation between simulta- 
neously active SSCs ^ 80 pc), a multiphase environment is 
built up in the central region of the galaxy favouring the 
formation of a filamentary-like wind, as in M82. 

The resulting filaments have extensions between ~ 50 
and 500 pc, thickness of ~ 30 pc, densities between 1 and 
10 cm -3 , temperatures of about 10 4 K, pressures between 
10~ n and 10" 12 dyne cm and velocities between ~ 100 
and 800 km/s. Outside of the disk, because the wind has 
a velocity ~2200 km/s, they are impacted by the hot gas 
with relative velocities of 1500-1800 km/s and therefore, the 
temperature may increase up to 10 7 K. Just a little fraction 
of these filaments escapes from the galaxy and, based on 
the results of this work, this corresponds to only ~ 10% of 
the gas in the high density phase, that is, about 5x10 M©. 
This result suggests that although the wind can be efficient 
to inject energy into the halo, it seems to be less efficient to 
inject disc mass (see also Section 4.4). 

We can conclude that the environmental conditions at 
the base will determine the main large-scale features of the 
GW morphology and particularly, whether filaments and a 
multi-phase environment will develop or not in the wind. If a 
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large number of SSCs form before the SNe driven hot super- 
bubbles fill out the entire nuclear volume of the galaxy, then 
an M82-like wind, with a rich filamentary structure, devel- 
ops. On the other hand, if only a few, highly massive SSCs 
form in the SB region then the SNe will drive a superbubble- 
like wind with no filaments and with all the phases of the 
gas (high, intermediate and low density) well separated from 
each other. 

The results above also suggest that the features of the 
wind mainly depend on the recent star formation history 
of the starburst rather than all earlier starburst episodes 
which occurred before ~ 50 Myr. In fact, although many 
observations indicate that intense episodes of star for- 
mati on occurred in M8 2 also about 100 Myr ago (see, 
e.g-. lde Griis et al]|200l| ), the simulations indicate that the 
present wind energy and morphology are essentially due to 
the recent activity occurring along few Myr of the younger 
SSCs lifetime. Remnants of fossil star formation activity had 
been swept long ago by the strong power driven by previ- 
ous SNe explosions and are possibly in the outer regions of 
the wind or have already evaporated. For this reason, it is 
sufficient to consider only the last ~10 Myrs of the star- 
burst lifetime in order to assess the main features of the 
large scale wind structure in the galaxy surrounds, as in the 
models discussed here. 



4.2 Energy budget 

The energy injected by the explosions of SNe from a sin- 
gle SSC has a heating efficiency HE ~ 100% (where, 
as remarked in Section 2.2, HE corresponds to the ra- 
tio between the kinetic plus thermal energies stored in 
the gas and the energy injected b y the SN explosions, 
iMelioli fc de Gouveia Dal Pinoll2004h . This is because the 
radiative losses are important only if t coo ; < t int (where t coo ; 
and ti n t are the gas radiative cooling time and the interac- 
tion time of the supernova remnants, respectively) and in the 
case of single SSC the SNe occur so close in space and time 
that this relation is almost never fulfilled. But what fraction 
of the total energy of the superbubbles goes to drive bulk 
gas motions and generate a large scale wind? In all the mod- 
els, we have found that the energy lost by radiative cooling 
is proportional to the ratio between the mass of the high- 
density gas component and the total mass of the GW. Since 
this value ranges between 70% and 82%, we verify that only 
a fraction between 20% and 30% of the SN energy goes into 
the kinetic energy of the gas to drive the GW (Figure [T3J. 
Therefore, in contrast to the results obtained for the for- 
mation of a multi-phase environment, the fraction of energy 
available for the build up of the GW is about the same for 
all the models independent on the temporal and space dis- 
tribution of the SSCs and on the SNe power injection during 
the first Myr. 

The injected energy that is not radiated away provides 
most of the kinetic energy to the gas in the high-density 
phase (filaments and clouds) and the thermal energy to the 
gas in the intermediate and low-density phases. These en- 
ergies are spatially distributed as depicted in Figure [T4l for 
model SC40-01. 

These results suggest that the average temperature of 
the hot gas will be given by: 



3 n i+ ik 

where e t h is the fraction of the total SNe energy stored in the 
gas as thermal energy and n i+ i is the total number density 
of the intermediate and low-density gas components. 

Similarly, the average velocity of the filaments will be 
given by: 

2e k E S N ,„-. 

Vf ~ \ (8) 

\ m h 

where e k is the fraction of the SNe total energy stored in 
the gas as kinetic energy, and mh is the total mass of the 
high-density gas. 

With the results of the Reference Model (SC40-01) pre- 
sented in Figure 1131 we obtain an average temperature for 
the low-density gas T ~ 10 7 K and an average velocity for 
the high-density gas Vf ~ 500 km s , where we have con- 
sidered e t h ~ 5% and tk ~ 18% (se e Figu re 20). This is 
in agreement with IChevalier fc Clegd (|l985l ) who predicted 
that for a thermalization efficiency of about 10% the wind 
fluid should have a temperature in a range between 10 7 
and 10 8 K even if it had been slightly mass loaded with 
cold ambient gas. We not e that although the wind model of 
IChevalier fc Clegd (|l985h is adiabatic, these authors intro- 
duced two parameters in their model (a and 0) in order to 
incorporate the effects of radiative cooling both from the SN 
ejecta and dense gas, so that a comparison of their model 
with t he present ones is appropriat e. However, as pointed 
out by Strickland & Hcckman (200^), such temperatures are 
hotter than the usually observed warm ionized, soft X-ray 
emitting plasma in GWs. Thus hard X-ray observations are 
still required in order to confirm or not the predictions of 
the models. 

Based on the results above of the energy distribution 
during the wind evolution, we can conclude that an M82- 
like GW could be generated in two steps. The first one can 
be characterized by the formation of supperbubbles due to 
SNe acti vity in each SSC with a SN heating efficiency HE 
~ 100% l|Melioli fc de Gouveia Dal Pindliooi '). The second 
step is characterized by large scale interactions of the su- 
perbubble shells with a HE (including the thermal and the 
kinetic energy contributions) of ~ 20% only. In this scenario, 
about 80% of the SN energy is radiated away mainly by the 
denser structures embedded in the GW. 



4.3 Chemical evolution 

Our simulations show that a large fraction of the metals in- 
jected into the ISM by the SNe are in the filaments located 
near the galactic disc (see the upper-right panel of Figure 
[BJ. In fact, the metal abundance in this high-density gas 
component is small, but the total amount of metals it con- 
tains is similar to the metal mass contents in the low-density 
component. 

Regardless of the model, the high-density, high- velocity 
structures have a maximum abundance close to the solar 
one, while the low-density phase is characterized by a maxi- 
mum abundance of ~ 10 Zq . Examining the maps of Figure 
[TT]for our Reference Model SC40-01, we can determine an 
average abundance, Z, for each gas phase. For the high- 
density phase we obtain Zh ~ 0.5 Zq, for the intermediate- 
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Figure 13. SN energy distribution for the different models stud- 
ied here after a time of 5 Myr. ft : fraction of the total SNe energy 
stored in the gas as vertical (or poloidal) kinetic energy; e^: frac- 
tion of the total SNe energy stored in the gas as thermal energy; 
e\: fraction of the total SNe energy lost by radiative cooling. 

density phase Zi ~ 1.2 Zq, while for the low-density phase 
we have Zi ~ 4.5 Zq. Each gas phase regulates the chemi- 
cal evolution of the different parts of the system. The high 
and intermediate-density phases, with their filaments and 
clumps, are unable to escape from the galaxy and are thus re- 
sponsible for retaining in it an important fraction of the met- 
als ejected by the SN explosions. In fact the total amount of 
ejected metals after 8 Myr is ~ 2 x 10 5 M Q , but the amount 
which escapes from our physical domain is only ~ 8x10 
Mq and the remaining portion is retained in and around 
the galactic disk. The total mass of metals in the high and 
intermediate-density phases are, respectively, Mh x Zh and 
Mi x Zi, that is, ~ 4xl0 4 Mq of metals for each phase, 
which represents ~ 2/3 of the metals stored in the galaxy. 
The remaining 1/3 of the chemical species is in the low- 
density phase and is therefore, the only fraction of metals 
that is transported to outside of the galaxy to enrich the 
IGM. 

Based on the results above, we can conclude that the 
SNe explosions in an M82-like SB galaxy change significantly 
the metallicity only of the low-density wind component. In 
our simulations this component has a resulting metal abun- 
dance between 5 and 10 Zq . The injected metals by the SNe 
do not affect significantly the metallicity of the high and in- 
termediate density phases of the wind, a result which is in 
contrast to what occ urs, for instance, in the cl ouds built up 
in galactic fountains jMelioli et al.ll200Sl , l2009h . 

4.4 Mass loss 

Our results suggest that during the early phases of expan- 
sion, the GW transports to above the disk a large amount of 
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Figure 14. Model SC40-01: map of the energy "column den- 
sity" (i.e., the energy density integrated in the transverse direc- 
tion to the wind) at a time 3.5 Myr. Top-left panel: total en- 
ergy; top-right panel: thermal energy; bottom-left panel: vertical 
or poloidal-kinetic energy; bottom-right panel: radial-kinetic en- 
ergy. The energy column density is expressed in erg cm~ 2 . 



energy, momentum and gas, but at the same time the high- 
density component (of clouds and filaments) is able to reach 
only intermediate altitudes and does not escape from the 
system. As a consequence, no significant amount of gas mass 
is lost in the IGM, except for an important fraction of met- 
als and SN ejecta that flows along with the high velocity-low 
density gas in the GW. In our Reference Model (SC40-01), 
for instance, the GW reaches a nearly steady state evolution 
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after 4-5 Myr and until that moment the mass lost by the 
system is about 0.6 Mq yr _1 only. In fact, as stressed in the 
previous Section (Sec. 4.3), the only gas phase which is able 
to flow freely into the IGM is the low-density phase. How- 
ever, due to its very low-density, the total gas mass in this 
phase represents only 5% of the total mass of the GW (see 
Figure [9]). We note th at this value is only sli g htly s maller 
than that indicated by IStrickland fe Heckmanl l|2009l ). 

Based on the results above, we can conclude that the 
mass evolution of the galaxy is not substantially affected by 
the SB events observed in the core of M82. 

5 SUMMARY 

The main conclusions of this work can be summarized as 
follows: 

(i) The environmental conditions at the base will deter- 
mine the main large-scale features of the GW morphology. 
If a sufficiently large number of super stellar clusters (SSCs) 
is build up before the SNe driven hot superbubbles fill out 
the entire nuclear volume of the galaxy, then an M82-like 
wind, with a rich filamentary structure, develops. In quan- 
titative terms, a multiphase wind is build up, as in M82, if 
the number of SSCs is > 7-10 while // is still // sC 1. 

(ii) About 80% of the SN energy is radiated away mainly 
by the high density structures which are formed from the 
fragmentation of the superbubble shells (see e.g., Model 
SC40-01). 

(iii) The SNe explosions in an M82-like galaxy change sig- 
nificantly the metallicity only of the low-density wind com- 
ponent. In our simulations this component has a resulting 
metal abundance between 5 and 10 Zq, while the metallicity 
of the high and intermediate density components is around 
the solar abundance (i.e., ~ 0.5 to 1 Zq). Also according 
to our results, about 1/3 of the total mass in metals is in 
the low-density component and therefore, this is the only 
fraction that can be removed to the IGM. 

(iv) During the early phases of expansion, the GW trans- 
ports to outside the disk large amounts of energy, momen- 
tum and gas, but the more massive high-density component 
(which is composed of clouds and filaments) is able to reach 
only intermediate altitudes (typically smaller than 1.5 kpc) 
and therefore, does not escape from the system. In conse- 
quence, no significant amounts of gas mass are lost to the 
IGM and the mass evolution of the galaxy is not much af- 
fected either by the starburst events occurring in the nuclear 
region. 
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